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Recent experiments in the cuprates have seen evidence of a transient superconducting state upon 
optical excitation polarized along the c-axis [R. Mankowsky et ah, Nature 516 , 71 (2014)]. Moti¬ 
vated by these experiments we propose an extension of the single-layer t — J — V model of cuprates 
to three dimensions in order to study the effects of inter-plane tunneling on the competition be¬ 
tween superconductivity and bond density wave order. We find that an optical pump can suppress 
the charge order and simultaneously enhance superconductivity, due to the inherent competition be¬ 
tween the two. We also provide an intuitive picture of the physical mechanism underlying this effect. 

Furthermore, based on a simple Floquet theory we estimate the magnitude of the enhancement. 


I. INTRODUCTION 

Ever since the discovery of high-temperature supercon¬ 
ductivity in cuprate material^, it has presented a funda¬ 
mental challenge to our understanding of strongly inter¬ 
acting systems. Despite the decades of strenuous efforts 
these materials still remain a profound mystery. In the 
last several year, however, new experimental results have 
shone some light on one of the most puzzling phases of 
cuprates - the pseudogapP. Clear evidence of some sort 
of charg e or dering competing with superconductivity has 
emerged^. 

Recently, in a series of exciting papers, it has been 
demonstrated that optical pumping can significantly 
enhance the superconductivity in several hole-doped 
cuprateP^^ (for a brief summary see Ref. [9]). In these 
works, a sample was subjected to a short pulse of mid- 
infrared light, polarized along the c-axis. Reflectivity 
measurements were then taken as a function of time de¬ 
lay, from which the frequency dependent conductivity 
was extracted. There is an increase in it that the au¬ 
thors identify with an inhomogeneous and transient en¬ 
hancement of superconductivity. Furthermore, the na¬ 
ture of this enhancement seems to be different from the 
already well known effect due to quasiparticle photo- 
excitatiorPESHm, Several works have investigated these 
experiments and proposed an incre ase of inter-layer cou¬ 
pling as one of the dominant effectHowever, there 
is another possible effect to consider, since resonant x- 
ray spectroscopy seems to suggest that the charge order 
is suppressed in the same region of the phase diagram 
where superconductivity is enhanced in the optical exci¬ 
tation experimentJi^. A natural question to ask, then, 
is how coexisting superconductivity and charge order be¬ 
have under such perturbation. 

To address this problem we start from the t — 
J — V model of the quasi two dimensional Cu02 


plane^i^ESI^ which can naturally support the coexis¬ 
tence and competition between charge ordering and 

superconductivitjESHSU Furthermore, we focus on the 

low energy physics of fermions near the so-called ‘hot 
spots’, where the F ermi surface intersects the magnetic 
Brillouin zon^^mS! y^e extend the model by introducing 
an effective Hamiltonian, describing stacked planes, cou¬ 
pled by a c-axis tunneling term Then we investigate 
the phase diagram of this extended model by utilizing a 
Landau expansion of the free energy. Quite surprisingly, 
we observe a non-monotonic behavior of the critical tem¬ 
peratures of the two orders with increasing tj., and we 
provide an intuitive physical explanation of this interest¬ 
ing feature. Finally, we consider different effects of the 
photo-excitations of the system, particularly focusing on 
the role of the apical oxygens , whic h are thought to play 

a key role in the experiment yyg that, quite 

generally, there is a parameter region where an optical 
pump can lead to a melting of the charge order and a 
corresponding enhancement in superconductivity, due to 
the competition between the two orders. 

As a side note, let us add that the t—J — V model 
and its variants tend to have a s their leading instabil¬ 
ity a (Q, Q) type ordering vectoiE^ISfllMI^ q) 

ordering vector seen in experimeniP as a sub-leading in¬ 
stability (with both orders having predominantly d-wave 
symmetry). While several extensions of the model have 
been proposed as a w ay to stabilize the experimentally 
observed ordeiEEHHSl^ they introduce additional, and for 
our purposes unnecessary, complications. In order to 
keep the model simple, we restrict our attention to the 
instability at wavevector {Q,Q). The physical content of 
our results lead us to expect that the qualitative behav¬ 
ior of the effects would be similar for the experimentally 
relevant {Q, 0) charge order. 

The structure of the paper is as follows. First, we pro¬ 
vide a short description of the two-dimensional t—J — V 
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model and the associated Landau theory in Section [TTj In 
Section mi we extend this model to include c-axis tun¬ 
neling, and explore the effects of on the phase dia¬ 
gram. Time-dependent perturbation of tz is introduced 
in Secti on |IV[ and the high-frequency limit is studied. In 
Section | V| we summarise and discuss our findings. 


II. TWO-DIMENSIONAL MODEL 

In order to study the interplay of bond density wave 
(BDW) and superconducting orders we e mploy a 2D 
t — J — V model of a Cu02 plan^l^lffSlMJ, it provides 
a natural platform for exploring the general features of 
the interaction and the coupling between these two or¬ 
ders within a single copper oxide plane. The Hamiltonian 
is 

H = JSi • Sj -l- - Vriirij, (1) 

(*j) (id) 

where V and J are nearest neighbor interactions, 
= So- is the charge density, and Si — 

^ the site spin density, with i and j being 

site indices, and cr, a, and b are spin indices. The term 
tij contains neare st ^ next to nearest, and next to next 
to nearest hoppin^^. V describes the nearest-neighbor 
tail of the Coulomb repulsion, which tends to suppress 
the d-wave superconductivity and enhance the BDW or¬ 
der. J is the usual nearest neighbor anti-ferromagnetic 
exchange interaction. 

Of course, the pure t — J model has been extensively 
used in the studies of cuprates as an effective one-band 
description of the Cu02 planed. It naturally leads to 
a d-wave superconductivity as its dominant instability. 
In order to have a region where superconductivity and 
BDW order coexist the model can be extended by the 
introduction of V, which suppresses superconductivity 
and boosts the charge order - this is the rationale behind 
the t — J — V model. 

We now construct a low-energy effective model by re¬ 
stricting our attention to fermions living within a limited 
region surrounding ‘hotspots’ - points where the Fermi 
surface intersects the magnetic Brillioun zone bound¬ 
ary, as depicted in Fig. (such models have been re¬ 
cently introduced and u sed in a n umber of studies of 
charge order in cuprate These points are 
of special interest because the Fermi surface is nested 
with wave-vector K = (7r,7r), and given the importance 
of anti- ferromag netic spin fluctuations to pairing in the 
cuprateJi^I^i^, we expect that the most relevant interac¬ 
tions will be those with exchanged momentum K. Close 
to the hotspots we can replace the interactions Vq with 
constants and restrict the two unconstrained 

fermion momenta k — k' = q to lie in hot regions sep¬ 
arated by K. 

Furthermore, we only consider a charge ordering in¬ 
stability with a diagonal {Q, Q) ordering vector and a 



FIG. 1. Hotspots exist where the Fermi surface intersects 
the magnetic brillioun zone boundary. The charge ordering 
vectors are given by the separations between hotspots 1 and 
2 accross the Brillouin zone border and the rotated vector 

{Q,-Q)- 


d-wave form factor, as this was found to be the strongest 
mean-field instabilitjff^HH jf further enforce time re¬ 
versal symmetry this allows us to concentrate our atten¬ 
tion to four hotspots. Having restricted our fermions to 
live within a range A of the hotspots we obtain an effec¬ 
tive Hamiltonian 
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the interaction is 


= - J^dab- Sad- ( 3 ) 


k is now the deviation from the hotspot, a — d are the 
electron spin indices, and i is now a hotspot index. 

At this point we undertake a mean field decomposi¬ 
tion of the interaction (four-fermion) terms of the Hamil¬ 
tonian simultaneously in the BDW and superconducting 
channels, where a d-wave form factor is assumed for both 
orders, by defining 


fcrL A = - 
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and 

(pSab = -gc =9 cY <4,fc.a^4,fc,b) ^ (5) 

k k 

where 

aAc4 

5c = 2 E 9^'’"‘'3adSbd = 3J + 4V, 

a,b,c,d 


are the effective couplings in the superconducting and 
charge channels respectively, and (■ ■ ■) indicates a ther¬ 
mal average. The order A is the projection of uni¬ 
form d-wave superconductivity onto the hotspots, and 
(p describes a d-wave charge order, that lives on the 
bonds between copper sites, and has modulation vec¬ 
tor Q given by the separation between hotspots. It is 
clear that V enhances superconductivity and simultane¬ 
ously suppresses charge order. Because of the d-wave 
form factor of the order parameters, two of the remain¬ 
ing hotspots become redundant and at the mean field 
level the behavior of the system may be described by a 
4x4 Hamiltonian in hotspot-Nambu space. Defining a 

Nambu spinor 
Hamiltonian takes the form 


H — 'tpiHMF{k)tpi: + 

- 9s 9c 

k 


where 


Hupik) = 
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( 8 ) 


and A and cp are the superconductivity and BDW order 
parameters respectively. 

Both order parameters are generally complex numbers, 
which we can write as A = |A|e*®* and (p = How¬ 

ever, since we can always remove the complex phases 
(at the mean-field level) via a gauge transformatiorP^ 
Cl —>■ C 2 —t we will con¬ 

sider only real and non-negative values for (p and A in 
our analysis. 

From [Equation 7| we can readily derive a Landau 
free energy for the A and (p orders. Evaluating / = 
fMF + ^ {H — Hmf)mf’ rising the above decoupling and 
expanding to fourth order in the order parameters, we 
obtain 


'2 

- ^BDW 

<P^ + 

2 

-Hsc 

V9c \ 


I9s J 


+ UBDW(p‘^ + USC^'^ + W(p'^A‘^, (9) 


with w,ubdw ,usc > 0. The exact expressions for the 
coefficients are given in Table I| It is in fact possible to 
solve the mean field problem exactly via a sequence of 
Bogoliubov transformations. This method breaks down, 
however, once we introduce c-axis hopping, and so, in 
anticipation of this extension of the model, we have cho¬ 
sen instead to work with a Landau expansion, which will 
carry over to the more complicated case. 


Term 

Expression 

Hsc 

TIbdw 

f l-2n(5i,T) 

Jk £i 

O r ri(£2,T)-n(£i,T) 

•Jk ^1—^2 

use 


ubdw 

fk [n'(Ci,T) +n'(6,T) + 2^5^] 

w 


TABLE I. Microscopic expressions for coefficients of the Lan- 


dau theory, n is the Fermi function and n' the derivative of 
the Fermi function with respect to energy. 
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FIG. 2. (Color online) Phase diagram of the 2D modepT!^ for 
fixed J. The regions are from left to right, SC, BDW -I- SC, 
and BDW, with the white region at the top being the normal 
phase. Note that the temperature does not extend down to 
T = 0 as the Landau theory is not valid in this limit. 

In what follows we hold J fixed and use V to adjust 
the splitting between the superconducting and BDW in¬ 
stabilities. Doing so, we obtain a phase diagram such 
as that in [Figure ^ In it we see the following easy-to- 
understand behavior: when V is zero or very small the 
d wave superconductivity is the only relevant instability 
of the model. With V increasing, superconductivity is 
suppressed and eventually BDW appears as the leading 
order. Notably, in the region where these instabilities are 
comparable there is a (rather narrow) coexistence phase. 
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These features are consistent with previous studies of the 
modeP^. 

We expect this phase diagram to be strictly valid only 
close to Tc and in the region in which the critical tem¬ 
peratures of both orders are comparable. However, com¬ 
paring the Landau expansion with the exact numerical 
mean-field solution shows that even for intermediate tem¬ 
peratures there is only a small, purely quantitative cor¬ 
rection to the shape of the coexistence region. 


III. EXTENSION TO STACKED PLANES 


Motivated by the experiments which have used optical 
excitation polarized along the c-axis to create transient 
states of enhanced superconductivitj^^, we seek to ex¬ 
tend the purely two-dimensional model from the previous 
section to include coupling between planes. The simplest 
model that encapsulates this behavior can be represented 
by the following Hamiltonian: 

( 10 ) 

I 

Here I is a layer index, and Hi is the single-layer Hamil¬ 
tonian considered in [Equation 1| There are copies of 
those, which are coupled via a c-axis tunneling term 

= E 4||.fc.>(^ii’*-)%.fe.,-’ (11) 

/C|| 


where c-axis momentum kz is conjugate to the plane in¬ 
dex 1. 

In what follows we will consider and compare three 
different forms for the tunneling tz in Eq. El The exact 
expressions for each tunneling type are presented in |Ta-| 


ble H Type A tunneling (nearest neighbor hopping along 


the c-axis) we introduce mainly for its simplicity, type B 
comes from a one-band ti ght binding fit to band struc¬ 
ture calculations for LSCCP^, and type C was proposed 
as an approximate tunneling form for several families of 
cuprate superconductoil^. Despite the significant differ¬ 
ences between these tunneling forms, it turns out that 
the effects we obtain are not specific to any of them, but 
are in all cases qualitatively similar. 


Type c-axis tunneling 

A —2tz cos kz 

—2tz cos (^) (cos kx — cos ky)^ cos (^) cos 
—2tz cos(kz){cos kx — cos ky)^ 

TABLE II. c-axis tunneling elements used in the calculation 


We now retrace the same steps as in [section H| The 
derivation proceeds similarly, but there are some sub¬ 
tleties that need to be considered first. 

Due to the model containing solely in-plane interac¬ 
tions, we only need to consider pairing of quasi-particles 


within the same plane. As a consequence of this, the or¬ 
der parameters do not depend on kz and the vector Q 
which separates fermions contributing to pairing in the 
charge channel cannot change with kz- Because of this 
restriction, we find that while at kz = 0 superconductiv¬ 
ity and BDW pair the same points in the 2D Brillioun 
zone, this is no longer true for kz ^ 0, tz ^ 0. The 
particle and hole being paired in the BDW channel must 
remain separated hy Q = (Q, Q) even when the Fermi 
surface is no longer nested with this vector, as can be 
seen in [Figure 3l Superconductivity on the other hand 


continues to pair k and —k for all kz- As a consequence, 
one has to be careful to define the hot regions properly 
in this case. There are several ways one might do so, but 
fortunately, as shown in the appendix, they all lead to 
the same qualitative behavior. 

That being the case, we implement the following pro¬ 
cedure. At each kz there is a region centered on where 
the 2D Fermi surface intersects the 2D magnetic Bril¬ 
lioun zone, (j) now pairs quasi-particles separated by the 
fixed charge ordering vector and is only non-zero when 
the momenta of both fall within a hot region. With this 
procedure in place, the free energy takes the same form 
as Equation 9 but with the coefficients now being the 


three dimensional integrals shown in [Table HI[ Here, we 
have made the simplifying assumptions that A and (j) are 
not modulated along the c-axis. As can readily be seen, 
the expressions in [Table HI[ reduce to those in [TableT] in 
the limit of no c-axis hopping. 



FIG. 3. c-axis hopping causes a warping of the Fermi surface 
which causes the hotspots to move as a function of kz and 
leads to a destruction of Q nesting away from kz = 0. 

Using each of the three tunneling forms, we calculated 
the state which minimized the free energy for a range 
of T, V and tz- The {Q,Q) order remains the leading 
charge instability at the quadratic level, so we again only 
decouple in this channel and the d-wave superconduct¬ 
ing channel. The phase diagram for tunneling type B is 
presented on Fig. (the other two types lead to quali¬ 
tatively similar diagrams). We start (for tz = 0) with 
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Term 

Hsc 

^BDW 

use 

ubdw 

w 




Expression 

r l-2n(gi7TJ 

Jk 

^ n(ei-Q,T)-n(gi,T) 

Sl-'Sl-Q 


[n'{^i,T) 


2/^ 


k 5l{ei-?l-Q) 


l-2ra(g 

25l 


iT) j 


TABLE III. Microscopic expressions for the coefficients of the Landau theory. ^i-q is the energy at a point at the point 
k' — k -\- kus — Q- Primed integration indicates a restriction of the integral to regions where both k and k' lie within a hot 
region. The kz integration is from —tv to n. 


the coexistence case in which BDW is the leading insta¬ 
bility. As we can see, for small but finite the charge 
order is further enhanced at the expense of superconduc¬ 
tivity. However, once becomes sufficiently large, the 
tendency reverses, and superconductivity is boosted by 
the increase of three-dimensionality, until it becomes the 
leading instability. 


4.0 



(meV) 


FIG. 4. (Color online) Phase diagram of a system of stacked 
planes at a fixed coupling strength for tunneling type B in |Ta-| 
|ble II| Types A and C exhibit qualitatively similar behavior, 
but with different scales on the tz axis. The inset indicates 
where the starting point at tz = 0 lies in E — T plane of 
Figure 2| {V = 110 meV). 


To understand this peculiar shape of the phase dia¬ 
gram, we look at the quadratic coefficients in the free en¬ 
ergy. The behavior of the superconducting coefficients, 
05(7 and uscj can be explained entirely by the depen¬ 
dence of the density of stat es and ba ndwidth, p{^, tz) and 
W{tz)- As can be seen in Figure 5 these lead to only a 
small effect on the superconducting susceptibility. Thus, 
the features visible in [Figure 4| are mostly determined by 
the behavior of the charge instability and its indirect ef¬ 


fect on superconductivity through the biquadratic term 
in the free energy. 

The non-monotonic behavior of the charge susceptibil¬ 
ity Ubdw can be understood as a competition of two 
effects. First, it is well known within this model that 
in the absence of Fermi surface curvature, the diagonal 
BDW order exhibits a BCS type transition characterized 
by a logarithmic divergence. Curvature however cuts off 
the diver gence of the logarithm and weakens the charge 
instabilitjffSESMl^ Specifically, we can write the disper¬ 
sions as 


Ci = c+ + ?- 

6 = C-i- - C- 


In that case the charge susceptibility becomes 


^BDW = 


f sinh ^ cosh 
k [sinh^ -I- cosh^ |^] 


( 12 ) 

(13) 


(14) 


in the 2D limit. One can see clearly that in the limit that 
—>■ 0, the absence of curvature, the charge (particle- 
hole) instability becomes as strong as the superconduct¬ 
ing (particle-particle) one, and for any non-zero the 
charge susceptibility is weaker compared to that of the 
BCS case. 

As the c-axis hopping tz increases the 2D curvature 
decreases. This is because the kz = 0 component for 

all tunneling types is of the form —tzTi{k^^)cl , with 

fc|i II 

rj > 0. Near the hotspots this acts as an effective upward 
shift in the chemical potential, which moves the hotspots 
such that is decreased relative to As a result, 
for increasing tz, the kz = 0 Fermi surfaces changes in 
such a way as to effectively enhance the charge instability. 
Fundamentally, this is a consequence of the effect of the 
tunneling term on the properties of the Fermi surface of 
each plane. 

In opposition to the aforementioned effect, increasing 
tz will lead to a progressive destruction of nesting away 
from kz = 0, as depicted in [Figure 3| For tz = 0, the 
Fermi surface is nested at the hotspots for all kz- How¬ 
ever, as tz increases the Fermi surface warps more with 
kz, decreasing the portion of the phase space for which 
there is approximate nesting, and thus weakening the 
charge susceptibility. 
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For small tz the warping of the Fermi surface is small, 
and so the 2D decrease of curvature is nearly the same for 
a wide range of kz values, leading to an overall strength¬ 
ening of the charge instability. However, as tz increases 
Fermi surface warping along the c-axis becomes more 
pronounced. As a consequence, the available phase space 
for charge ordering is substantially reduced, while at the 
same time the significance of the decreased curvature is 
lessened away from kz = 0, together eventually caus¬ 
ing a weakening of the charge ordering instability. The 
non-monotonic shape of the phase diagram can conse¬ 
quently be understood as demonstrating a crossover be¬ 
tween regimes in which the 2D and the 3D effects of tz 
dominate. 



enhancements of tz will generally lead to suppression of 
charge order and enhancement of superconductivity. 

It should be noted that the effects demonstrated here 
are obtain within the grand canonical ensemble at fixed 
chemical potential. In anticipation of section fV| the ex¬ 
perimental program we envision comprises placing a sam¬ 
ple in contact with its environment, which functions as 
a particle bath maintaining /i, and then applying pertur¬ 
bations that will lead to a change of tz in the effective 
Hamiltonian. This, however, means that the volume en¬ 
closed by the Fermi surface is not constant and therefore 
average particle number is not conserved. If one were in¬ 
stead to consider a system where such a constraint were 
important this might change the behavior of the phase 
diagram at small tj., where the effects are largely gov¬ 
erned by the position of the kz = 0 Fermi surface. How¬ 
ever, larger tz should still lead to Fermi surface warping 
and thus the enhancement of superconductivity in that 
regime should remain. Therefore, while the 2D effects 
that are important predominantly at small tz could be 
washed away, the suppression of charge order (via de¬ 
struction of nesting), and thus the enhancement of su¬ 
perconductivity, appears much more robust. 


IV. ENHANCEMENT OF 
SUPERCONDUCTIVITY VIA PERIODIC 
MODULATION 


(a) Quadratic Susceptibilities 





0 50 

u (meV) 


(b) Quartic Terms 

FIG. 5. (Color online) Quadratic susceptibilities and quar¬ 
tic coefficients of the Landau free energy as a function of tz, 
scaled by their value at tz = 0, for type B tunneling (at fixed 
temperature). Types A and C again exhibit similar behavior. 
Note that there is no significance to the crossing of lines for 
different coefficients, as the leading instability will be deter¬ 
mined by ^n. 


Interestingly, the value of tz for which superconduc¬ 
tivity reaches a minimum is of similar magnitude to the 
strength of c-axis hopping for various families of cuprates 
obtained as a fit to band structur^^. Given this, if we 
imagine a system where tz is near or greater than the 
point of minimum superconducting Tc in Figure 4 then 


Very recently, it has been argued that in YBCO the c- 
axis vibrations, induced by the external field, can change 
the equilibrium lattice structures). The result is a tran¬ 
sient shift of the Cu 02 planes - the intra-bilayer distance 
increases, while the inter-bilayer one decreases. In our 
model this would lead to effectively bringing the layers 
closer. Intuitively, it is clear that this should lead to en¬ 
hancement of the inter-plane coupling tz ■ In this the role 
of the apical oxygen seems quite significant; previous ex¬ 
perimental works have found that there is a range of dop¬ 
ings for which the inter-layer hopping exhibits a roughly 
exponential dependence on the hole-doping^ while the 
bond length between the plane copper and the apical 
oxygen exhibits a roughly linear dependence on doping 
through the same regioiP^. A natural interpretation is 
that the hopping exhibits an exponential dependence on 
an effective barrier width d\ 

U = Ae-“^ (15) 

where d is approximately linear in th e distance between 
the plane and the apical oxygerP^l^s) Thus, decreasing 
the interlayer distance leads to an exponential increase of 
the c-axis tunneling. 

There is a second, more subtle, way in which driving 
the apical oxygens can enhance the interlayer tunneling. 
Let us consider a vibration of these ions without change 
of their equilibrium positions. Then, in line with the 
above reasoning, we can model the effect of harmonic 
oscillations of the apical oxygen on tz as oscillations of 
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the effective barrier width, leading to a time dependent V. DISCUSSION AND CONCLUSION 

hopping element 


= Ae exp[—adi cos(nt)], (16) 


where d{t) = do + di cos(Ot), and tzo = Ae 

Let us consider the high frequency limit. Then, we 
expect that the quasiparticles will see an effective time 
averaged Hamiltonian. While in experiment the frequen¬ 
cies are not extremely high, we consider this limit as 
a particularly simple case, from which we may extract 
relevant qualitative trends. More formally, we can ob¬ 
tain a Floquet Hamiltonian related to the time depen¬ 
dent hoppings. The stroboscopic dynamics of the system 
will be governed by this Floquet Hamiltonian, which can 
be obtained as series in 1/H via a Magnus expansion of 
the time evolution operatoif^. For a high frequency os¬ 
cillation, we keep only the first term of this expansion, 
which is just the time-dependent Hamiltonian averaged 
over one period. In this case the Floquet Hamiltonian 
is the original Hamiltonian (Equation 101, but with the 
modification 


f ^ dt 

tz -t {tz{t)) = n J = tzoIo{adi), (17) 

where /q is the modified Bessel function of the first kind. 
lo is bounded below by 1, and increases monotonically 
with the magnitude of c?i. Therefore, within this approx¬ 
imation, any oscillation of the apical oxygens unavoidably 
leads to enhancement of the effective c-axis tunneling tz- 
This can be easily understood from the exponential de¬ 
pendence of the tunneling on the apical oxygen position: 
in the tail of the exponent, a stronger enhancement is ob¬ 
tained from decreasing the argument than the suppres¬ 
sion when increasing it at tt/VL time later. Thus, the 
tunneling amplitude is, on average, enhanced. Note that 
this observation is rather general, and likely applicable 
well beyond the region of validity of the high-frequency 
approximation used above. 

By using data from Refs.[37land|38]we can obtain an es¬ 
timate to the magnitude of both the oscillatory enhance¬ 
ment and the slower, quasi-static, effecl!^. Assuming an 
oscillation distance of di ^ 2.2 pm, as is observed in the 
experiments discussed in Ref. M we find approximately 
a 10 — 15% enhancement of tz in the steady state aver¬ 
age. This in turn may lead to a few percent up to around 
a 15% enhancement in the superconducting Tc depend¬ 
ing where in the non-monotonic structure of Fig. the 
sample is before perturbation. A stronger effect is the 
shift of the equilibrium position of the apical oxygen. A 
quasi-static shift of the apical oxygen position by 2.4 pm 
(again see Ref. [14]) can enhance Tc by up to about 60%. 
These effects provide a qualitative picture of a possible 
mechanism underlying the explanation given in Ref. 1141 


The primary motivation for this investigation came 
from the recent experiments on transient enhancement 
of superconductiv ity in the cuprates via mid-infrared 
optical excitation^SHUm Xo model these experiments 
we considered an extension of the t — J — V model of 
cuprates to three dimensions, and the effects of this three- 
dimensionality on the competition between superconduc¬ 
tivity and bond density wave orders. We showed that 
for our extended model, increased inter-plane tunneling 
leads to a suppression of charge ordering, and a coincid¬ 
ing enhancement of superconductivity due to the inher¬ 
ent competition between the two orders. The evolution 
of charge order takes place in two steps, corresponding 
to the regions where 2D and 3D effects of increased in¬ 
terlayer coupling, respectively dominate. The primary 
effect of interest is that the charge instability is sensi¬ 
tive to the c-axis curvature of the Fermi surface, which 
destroys nesting at the charge ordering vector. This ef¬ 
fect is generic across several tunneling forms proposed for 
various cuprate materials. 

These results provide a physical picture explaining the 
enhancement of superconductivity by the decrease of the 
inter-bilayer distance, caused by optical excitation. We 
further showed that periodic oscillations of the apical 
oxygens (identified in the experiments as important) can 
also lead to an effective increase of the inter-layer cou¬ 
pling. Both these effects indirectly promote superconduc¬ 
tivity via a suppression of the competing charge ordering. 
We believe that the mechanism presented here could play 
a significant part in the observed enhancement of super¬ 
conductivity and could be useful in pursuing new ways 
to raise Tc- Other mechanisms have been proposed with 
regard to t hese experiments: the suppression of phase 
fluctuationJi^EIl .^^gll as the usual enhancement of su¬ 
perconductivity due to microwave stimulatioiPSl could 
certainly play a complementary role. 

At the end, let us also note that this work considers 
only the mean field behavior of such a system. As is well 
known, fluctuations play an extremely important role in 
the superconducting transition of cuprate^^SHMl Never¬ 
theless, we expect that the effects discussed here on a 
mean field level will remain important in a more com¬ 
plete description of the system (entering through the rel¬ 
evant energy scales, for example). The mean field theory 
presented in this paper is only the first necessary step in 
the study of these effects, and including fluctuations is 
an important direction for future work. 
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Appendix: Low energy model away from kz = 0 



grations over the same hot regions as the superconducting 
terms, but enforce that pairing occur between particles 
separated by the fixed ordering vector Q. We label this 
the moving approximation. 

A final approach, and the one we used in this work, is 
to take a similar tact as the moving approximation, but 
to further restrict the integration such that both of the 
paired particles lie in a hot region. We will call this the 
truncated approximation. If we are committed to the no¬ 
tion that only fermions near the hotspots are important, 
this seems the most natural of the approximations, as we 

— Unmoving 
- - Moving 




FIG. 6. Different choices of hot regions for charge integrals. 
The gray dashed lines denotes the unmoving hot regions. The 
black dashed lines represent the moving hot regions, where 
one must be careful to avoid double counting. Finally the 
gray-hlled areas represent the truncated hot regions. 

As discussed in Isection IIII there are subtleties associ¬ 
ated with defining the hotspot model for stacked planes 
with inter-plane hopping. Specifically, one has to define 
the integrals involving charge order such that the separa¬ 
tion between paired particles remains constant with kz, 
despite the fact that the Fermi surface changes shape. 
We considered three different methods of handling this 
issue. Each is depicted in [Figure"^ In all cases the coef¬ 
ficients 115(7 and Use are unchanged, so we only need to 
decide how to implement the other three: IVbdw, ubdw 
and w. 

The first is to define two different hot regions: one re¬ 
mains bound to the Fermi surface and is associated with 
superconductivity, the other corresponds to the kz = 0 
hot regions for all kz and is associated with charge order. 
All integrals involving charge order are done over the un¬ 
moving hotspots. We label this procedure the unmoving 
approximation. 

Another approach is to perform the charge order inte- 


0 50 

(meV) 


FIG. 7. (Golor online) Quadratic susceptibilities as a function 
of c-axis hopping for different choices of integration region in 
integrals involving charge order. The qualitative behavior is 
the same for all three schemes. 


are then only counting the fermions that live within the 
true hot regions. The coefficients ^bdw, ubdw, and 
w are shown for a range of tz in [Figure 7 As can be 
seen, while the results differ numerically, the qualitative 
behavior is the same. Since the hotspot model is itself a 
qualitative model, we chose to adopt the truncated ap¬ 
proximation as it seemed the most in line with the spirit 
of the model. 

While these issues make the 3D model slightly more 
complicated, the qualitative behavior of the system ap¬ 
pears to be fairly robust to their method of resolution. 
This suggests that the extended model, like its simpler 
progenitor, can be used as a tool to uncover general phys¬ 
ical mechanisms underlying the complicated behavior of 
various lattice models. 
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